Signature of Quantum Hall Effect Skyrmions in Tunneling: A Theoretical Study 
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We present a theoretical study of the I — V tunneling characteristic between two parallel two- 
dimensional electron gases in a perpendicular magnetic field when both are near filling factor 
v — 1. Finite-size calculations of the single-layer spectral functions in the spherical geometry and 
analytical expressions for the disk geometry in the thermodynamic limit show that the current in 
the presence of skyrmions reflects in a direct way their underlying structure. It is also shown that 
fingerprints of the electron-electron interaction pseudopotentials are present in such a current. 

PACS numbers: 73.23.Hk 
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Over the years, tunneling experiments have proven to 
be a powerful source of information in many-body sys- 
tems. The two-dimensional electron gas (2DEG) in the 
quantum Hall regime |l|] is no exception. Both edge and 
bulk tunneling densities of states have been probed in 
different regimes of the bulk filling factor z/ = N/N^ [N 
is the number of particles and N^ is the number of mag- 
netic flux quanta). Recent tunneling measurements [0 on 
edge states have shown, for example, the first indications 
of chiral Luttinger liquid behavior, which has been the 
focus of much theoretical |^ and numerical Q study. For 
bulk systems, a correlation pseudogap observed in tun- 
neling experiments H has also motivated a large body of 
theoretical work ||] . 

The region very near v — \, however, has received 
much attention recently since the quasiparticles emerging 
from adding or removing charges turn out to carry spin 
textures, and are commonly called skyrmions [0-0]. In a 
simple characterization they consist of a spin-polarized 
quasiparticle with K (an integer) additional electrons 
flipped into the minority direction ||l3-n3]. This in turn 
may be interpreted as a charged quasiparticle with K 
spin- waves bound to it. The precise value of K depends 
on the details of the electron-electron interaction in the 
2DEG as well as the strength of the Zeeman coupling. 
It thus may be varied by tilting the sample with respect 
to the direction of the magnetic field [|l^, or by more 
involved methods [|5|. A strong experimental case that 
excess spins arc flipped when the quantum Hall ferro- 
magnet is doped away from v = \ has recently emerged 
|l4| , [l5[ , with the corresponding value of K in remarkably 
good agreement with theory |Q,|lJ,|l^. However, these ex- 
periments do not actually show that the flipped spins are 
bound to the excess charges. In this work, we show that 
tunneling experiments can directly demonstrate this ba- 
sic property of skyrmions, and measure, in a simple way, 
both the binding energy and number of spin waves bound 
to the quasiparticle. 

We consider a system in which there are two parallel. 



weakly coupled, 2DEG's for which there are dilute and 
equal densities of quasiparticles around the v = \ state. 
For the low densities considered, a treatment in which 
inter-skyrmion coupling |lj] is neglected is appropriate. 
Coulomb repulsions tend to keep the skyrmions in the 
same layer and opposite layers apart, so that the problem 
reduces to considering tunneling between a skyrmion and 
the V — 1 ferromagnetic state. Disorder effects are also 
ignored in our calculations. 

In addition to energy conservation, the conservation of 
spin greatly limits the possible tunneling events that can 
occur. For example, no tunneling is possible between 
regions of pure v = 1 ferromagnet because no unoc- 
cupied states are available for the lowest Landau level 
(LLL), spin majority electrons to tunnel into B. Tun- 
neling can only occur in the vicinity of the quasiparti- 
cles. (In what follows, we will consider for concreteness 
only quasihole excitations, i.e., v < 1. The results for 
u > 1 are identical due to an exact particle-hole symme- 
try pl p^jl^ .) For the case of spin-polarized quasiparti- 
cles {K = skyrmions), this can only occur precisely at 
zero voltage {V = 0) between 2DEG's. Thus, one ex- 
pects a 5 function response in the I ~ V curve which is 
presumably broadened into a narrow peak by weak dis- 
order. 

The case oi K = 1 quasiparticles, which is the major 
focus of this work, is richer and far more interesting. Fig. 
|l| illustrates the I —V characteristic as calculated in sev- 
eral different ways, and summarizes the central results of 
this work. The dashed line corresponds to an exact result 
obtained from a finite-size calculation of electrons on a 
sphere. A small broadening has been introduced into the 
(5-functions (which one always gets in any finite-size cal- 
culation) to mimic the expected continuous curve of the 
thermodynamic limit. In sharp contrast to the K = Q 
case, the tunneling current has a peak at finite V , and is 
suppressed near V = Q. However, in order to understand 
the nature and functional form of this supression, it is 
necessary to take the thermodynamic limit. As will be 



shown below, this can be done analytically in the disk 
geometry. The result obtained from using a simple wave- 
function for the K = 1 skyrmion |1J] is shown as a solid 
line in Fig. y. The tunneling current is completely sup- 
pressed below a threshold voltage, and then jumps to a 
finite value. This jump, however, turns out to arise from 
the fact that such a wavefunction is not an eigenstate of 
5^ . We propose below a new K = 1 wavefunction that is 
an eigenfunction of 5^ . The resulting I — V curve rises 
now linearly from the threshold voltage as illustrated by 
the dotted line in the lower inset. 

To understand the origin of the threshold, it is useful 
to consider the forms of the initial and final states in the 
tunneling process [see panel (a) in the upper inset]. A 
majority spin electron tunnels from the ferromagnet and 
fills one of the holes of the K = 1 skyrmion, liberating 
one particle-hole pair. Note that the hole left behind 
in the ferromagnet is completely equivalent to a i^ = 
skyrmion state, i.e., a spin-polarized quasiparticle. The 
lowest energy state that can be a final state on the right 
side is the fc — > spin- wave, which has an interaction 
energy that is identical to the ferromagnetic state. Thus, 
the minimum energy that one must supply is the differ- 
ence in interaction energy between the K — Q and K = 1 
skyrmion states Uk=o — Vk=\- This is the threshold en- 
ergy cVt- 

We will argue below from symmetry considerations 
that for general values of iiT > there will be a gap in the 
/ — F, and that the onset of current above the threshold 
takes the form !{¥) ex [F - Vrf^-^ for V >Vt. Here 
Vt is the interaction energy gained by deforming a spin- 
polarized quasiparticle into a skyrmion of size K [n3J16| 
{Uk=o—Uk)- Thus, tunneling spectroscopy measures ba- 
sic properties of the skyrmion in a direct way: the voltage 
threshold for the onset of current is the binding energy 
of the spin texture to the quasiparticle; the power law 
with which the current first begins to flow measures the 
spin of the skyrmion. For experiments on non-disordered 
samples, this implies that the tunneling I — V charac- 
teristic changes sharply as the Zeeman coupling is varied 
(as for example in tilted field geometries), whenever the 
equilibrium value of K changes. 

Finally, the high-energy behavior of the tunneling /— V 
curve displays sharp peaks for K = 1 skyrmions (Fig. 
n[). These arise from tunneling processes in which the 
minority electron tunnels out of the skyrmion into the 
ferromagnet [panel (b) in upper inset]. The final states 
on the skyrmion side are now polarized quasihole pairs, 
whose eigenenergies are determined only by the relative 
angular momentum of the two holes and are given by 
the pseudopotential parameters Vj M of the 2DEG. It 
is worth noting that if these peaks can be observed, this 
would constitute a direct measurement of Vj. To our 
knowledge, there are no other such direct methods for 
measuring these basic interaction parameters. 

The tunneling current between the two parallel 



2DEG's is given to lowest order in perturbation theory 
by 

/oo 
dujA-,AoJ - eV/2)At,{u + eV/2), (1) 

where the sum runs over all the single particle states in 
the LLL. (Higher Landau levels are not considered in this 
work.) Ama are the single-layer spectral functions for 
adding (removing) a particle to (from) a state m with 
spin a. In order to obtain the contribution to the cur- 
rent coming from the type of processes depicted in Fig. 
1, panel (a), we need four spectral functions for each m: 
^mT'^^T'^mi' ^^^ ^™i- Taking t as majority spin, 
A^^, = 0, which leaves us only with the spin-majority 
processes ^^T^mf ^'^^ ^^^ ferromagnetic state and an 
appropriately chosen zero of energy, A~^Juj) = S{uj) and 
we only need to calculate A'^^Jlj). To get the qualita- 
tive behavior of this spectral function we first perform 
finite-size calculations on a sphere [Q]. For each value of 
the total spin operator S^ , single skyrmions are obtained 
when N = N^ - I ^d^]. We will focus on the sim- 
plest one, S — N/2 — 1 {K = 1). When such a skyrmion 
receives a spin-majority electron the resulting state is a 
spin wave whose form is known exactly in the spherical 
geometry |l^], so that calculation of A^Jco) is straight- 
forward. The result of this calculation for A^ = 30 corre- 
sponds to the dashed line shown in Fig. |l|. As expected, 
the major contribution comes from injecting the electron 
at the center of the skyrmion which gives rise to the peak 
around V = OAe^/£. 

Because A'^Juj) is small for small cu, the behavior of 
the / — V^ at low voltages cannot be distinguished from 
finite-size studies. In order to address this we turn our 
attention to the disk geometry. The spin waves states of 
momentum k are given by 
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where the c^»,|s's create electrons in guiding center co- 
ordinate states, and |^o) represents the ferromagnetic 
state. The magnetic length [i = ^/hc/JeB)] serves as 
our unit of length. The corresponding energies, e(fc), are 
also easily evaluated |1^]. We need to compute the over- 
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lap of these states with cl^JK = 1), where \K = 1) 
represents the skyrmion state. In general, the neutral 
object Cm^K — 1) may always be written in the form 
I^me.m,. "™emhl'TiT;"te"ift), i.e., as a sum of states with 
one hole in the majority spin (nih) and one electron (we) 
in the minority spin. Written in this form, the matrix el- 
ement (fc|cjjj||-ftr = 1) may be computed exactly, and, 
furthermore, the sum over all final spin wave states in 
the spectral function can be computed analytically in the 
thermodynamic limit. The result of this calculation takes 
the form 



A+ H = fcoe-'=o/^ 
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where fco is impUcitly given by e(fco) = w, P'^^^ra2 

'mi\ f 1712^ 
n J \ n 
largest integer for which l3^-^„^^ is defined. 

In order to use Eq. ||, one needs to have an exphcit 
form for the skyrmion wavefunction, which supphes the 
coefBcients {a}. For the K = 1 skyrmion, the simplest 
variational wave function that has -?C as a good quantum 
number looks schematically like |12] 
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where the upper (lower) circles represent the spin- 
minority (-majority) single-particle states in the circular 
gauge, and the optimal coefficients have been found to 
be given by 7^ oc " "^'' |l|]. 

The values of amlnih can be obtained straightforwardly 
for this wavefunction, and the resulting I—V is presented 
in Fig. p^ (solid line) . It is interesting to note that in this 
approximation we find A^Juj) oc w'™"^', so that all the 
spectral functions vanish as w — » except for the case 
m — 1, leading to the jump in current at Vt- How- 
ever, a symmetry analysis shows that this behavior can- 
not be correct: all the spectral functions A^Jui) must 
vanish as w -^ 0. To see this, consider the matrix element 
{k\cl^JK = 1). The spin wave operator S~{k) (Eq. ||) 

commutes with cl^^ so that this expression may formally 

be rewritten as {K = 0|S'"(fc)|A' = 1) {\K = 0) repre- 
sents a polarized quasihole state). In the limit a; ^ 
only the k —> matrix element enters A^^{uj). How- 
ever, S~{0) is just the spin lowering operator, which can 
change only the Sz eigenvalue of any eigenstate of the 
Hamiltonian, but not its total spin S eigenvalue. Since 
skyrmions of different K values are in different 5^ multi- 
plets 0, it follows that {K = 0\S~{k)\K = 1) vanishes 
as fc — > 0, presumably as fc, the lowest order allowed by 
symmetry. Noting that the spin wave spectrum of the 
ferromagnet has energy state oc k^ for small k, one can 
easily see from Eq. that all the spectral functions must 
vanish at least linearly with w for the K = 1 skyrmion. 



That Eq. ^ above does not respect this property is due 
to the fact that it is not an eigenfunction of S"^. We pro- 
pose an improved wavefunction that is in the appropriate 
S"^ multiplet: 
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with the condition 71 + X]m=i I'm — guaranteeing that 
Eq. is in a well-defined S"^ multiplet. The coefficients 
{7m J 7m} can be obtained numerically by diagonalizing 
the Hamiltonian in such a Hilbert subspace. Table I il- 
lustrates the overlaps of the two approximate wavefunc- 
tions with exact-diagonalization wavefunctions on a finite 
disk. Although both states have very large overlaps, the 
\K = I)-)., 7' wavefunctions are clearly two orders of mag- 
nitude closer to the exact wavefunctions. The dotted line 
in lower inset shows the I-V curve at the onset of current 
for \K = l)-^,^' with the coefficients extrapolated to the 
thermodymanic limit. The current vanishes linearly near 
Vt since now A^^ has the proper behavior; the rest of the 
curve is essentially identical to the result obtained using 
the \K = l)j wavefunction. 

The contribution to the I-V characteristic coming from 
the spin- minority processes shown in panel (b), upper 
inset, dominates the voltage region above ^■njle'' jl. 
When a minority-spin electron is removed from a A" = 1 
skyrmion two quasiholes are left behind. The eigenstates 
and energies of this system have well-defined relative 
(J=odd integer) and center of mass angular momentum; 
the energy of the state depends only on J and is given 
by the Haldane pseudopotential parameters Vj 0]. The 
weights associated with these states are easily calculated; 
for the lAT = 1)^ one finds 
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The results are shown in Fig. H^ where the delta functions 
have been broadened to make them visible. 

What about skyrmions with larger values of AT? Ex- 
act analytic expressions for the spectral functions can- 
not be obtained using the methods above. However, we 
can draw some conclusions about what they and the re- 
sulting I — V must look like. Based on the previous 
reasoning, there will be a gap in the tunneling I — V, 
with current starting to flow at eVr = Uk=o — Uk- For 
voltages just above Vt, the final state on the skyrmion 
side may be interpreted as a linear combination of K 
spin waves, all with small wavevectors. Such states 
may be written in the form rii=i '^"(^Ol^o)- Noting 



again that c!^^ commutes with S (fci), the matrix ele- 
ments entering the spectral functions may be written as 
(K = 0| n^i S-{ki)\K). All the fci's must be non-zero 
in order for this expression to be non- vanishing, just as 
in the case described above for if = 1. Thus, the matrix 
elements entering into Eq. ^ should be proportional to 
Yl- kf . Assuming that the spin waves may be treated as 
non-interacting in the long-wavelength limit, the energy 
of the K spin wave state will be oc ^^ kf. Using these 
observations, one finds that A^Jcj) must vanish at least 
as fast as (w— Vt)'^^~^- The resulting current then rises 
from Vt as {V — VtY^^^ ■ Thus, for any skyrmion with 
K > 1, the threshold voltage and power law for the onset 
of current are respectively measures of the binding energy 
of the spin texture and the total number of flipped spins 
bound to it. 
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FIG. 1. Tunneling I-V in the presence of skyrmions for 
disk (solid line) and spherical (dashed line) geometries (see 
text.) Upper inset: Schematic representation of the two pos- 
sible tunneling processes, (a) A majority-spin electron tunnels 
into a skyrmion. (b) A minority-spin electron tunnels from 
the skyrmion into a ferromagnetic region. (Only the domi- 
nant Slater determinant is shown.) Lower inset: Blow-up of 
the low voltage region of / — V using simple (solid line) and 
improved (dotted line) wavefunctions. See text. 



